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ABSTRACT 

The 11 known satellite galaxies within 250 kpc of the Milky Way lie close to a great 
circle on the sky. We use high resolution N-body simulations of galactic dark matter 
halos to test if this remarkable property can be understood within the context of the 
cold dark matter cosmology. We construct halo merger trees from the simulations and 
use a semianalytic model to follow the formation of satellite galaxies. We find that 
in all 6 of our simulations, the 11 brightest satellites are indeed distributed along 
thin, disk-like structures analogous to that traced by the Milky Way’s satellites. This 
is in sharp contrast to the overall distributions of dark matter in the halo and of 
subhalos within it which, although triaxial, are not highly aspherical. We find that 
the spatial distribution of satellites is signihcantly different from that of the most 
massive subhalos but is similar to that of the subset of subhalos that had the most 
massive progenitors at earlier times. The elongated disk-like structure delineated by 
the satellites has its long axis aligned with the major axis of the dark matter halo. We 
interpret our results as reflecting the preferential infall of satellites along the spines of 
a few filaments of the cosmic web. 


1 INTRODUCTION 

In the cold dark matter (CDM) cosmology, structure builds 
up through fragments merging together in a roughly hierar¬ 
chical way. High resolution N-body simulations of the for¬ 
mation of dark matter halos in the ACDM cosmology have 
demonstrated that the cores of tightly bound fragments of¬ 
ten survive the merging process and re main as distinct sub¬ 
structures orbiting inside a parent halo iKlvoin et a Dim 
iMoore et alJI 19931 . The centre of the main halo and the ac¬ 
companying substructures are naturally identified with the 
formation sites of central and satellite galaxies respectively. 
The N-body simulations suggest that the mass functions 
of surviving substructures in galactic and cluster halos are 
roughly self-similar. Yet, the luminosity function of galaxies 
in rich clusters has a very different shape from the luminos¬ 
ity function of satellites in s maller systems such as the Milky 
Way or the Local G roup Jkauffmann et al.l Il993l: i Mated 

ll998l : lTrentham fc Hodgkinl200^lBenso^^^inl2nn2ai ~Not 

only is the central galaxy much more prominent in Milky- 
Way systems than in galaxy clusters, but the number of 
surviving subhalos in simulations of galaxy-sized halos far 
exceeds the number of the known satellites of the Milky 
Way. 

The discrepancy between the small number of satellites 
around the Milky Way and the large number of surviving 
substructures, once regarded as a major challenge to the cold 
dark matter cosmology, is now thought to be due to the as- 
trophysical processes that regulate the cooling of gas in halos 


and its subsequent transformation into stars. The increase 
in the entropy of the intergalactic medium brought about by 
the reionization of the gas at early times has been identified 
as a possible solution to the so-called “satellite problem” 
iKauffrnmin_ei^_alj|lT9jjBulIockl l2002t iBenson et al.ll2002al : 

see also lstoeh^^dl^OoJi . R.eionization sharply reduces the 

efficiency of gas cooling in small halos so that galaxies that 
formed prior to reionization are preferentially those that end 
up as satellites in syst ems like the Local Gro up. The detailed 
model calculated by iBenson et a 0 (l2002a^ which includes 
the effects of early reionization as well as other forms of 
feedback, reproduces many observed properties of the Local 
Group’s satellite system, including the distribution function 
of circular velocity, the luminosity function and the colour 
distribution. 

While the original satellite problem is no longer deemed 
a serious challenge, another related potential problem for 
the cosmological paradigm has recently been highlighted by 
Kroupa, Thies & Boily (2005). These authors argue that 
the strongly flattened spatial distribution of the 11 brightest 
dwarf satellites of the Milky Way, a feature known, but not 
unde rstood, for many years llLvnden-BelJ llOsD : iMaiewsl^ 
[l99^, is inconsi s tent w ith the ACDM model. According 
to iKroupa et all (( 200 D), CDM models predict a roughly 
isotropic distribution of satellites. They based this con¬ 
clusion on the assumption that the spatial distribution of 
satellites resembles the spatial distribution of the halo dark 
matter which indeed, as N-body simulations have demon- 
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strat ed, is approximately (although not exactly) spher ical 
(e.g. iFrenk et, alJIlQSst I.Ting fc Sutoll20o3 : lBulloclJl2002^ . 

In this paper, we demonstrate that the satellites of sys¬ 
tems like the Local Group do not trace the distribution of 
halo mass. On the contrary, the satellites in our suite of 
high resolution N-body simulations are generally arranged 
in highly flattened configurations which have similar prop¬ 
erties to those of the Milky Way satellite system. This, at 
first sight surprising, result is a reflection of the anisotropic 
accrection of subhalos which generally stream into the main 
halo along the filaments of the cosmic web. The flattened 
structure in which the Milky Way satellites lie traces a great 
circle on the sky and is approximately perpendicular to the 
Galactic Plane. In our simulations, the satellite structures 
tend to be aligned with the major axis of the triaxial halo 
mass distribution, that is, the longest axis of the halo is close 
to lying in the principal plane of the satellite distribution. 

As this paper was nearing completion, two related 
papers appeared on astro-ph. Both of them used high- 
resolution simulati ons of galaxy halos similar to those that 
we have performed. iKane et alJ identified “satellites” 

in their 4 simulations with randomly chosen dark matter 
particles taken either from the halo as a whole or exclu¬ 
sively from substructures. They were able to hnd flattened 
satellite systems similar to that of the Milky Way in t he for¬ 
mer case but not in the latter. IZentner et alJ ll2^^!^^ found 
satellites in three N-body simulations of Milky-Way type 
systems also in two different ways. In the first, they used the 
semi-analytic model of Kravtsov, Gnedin & Klypin (2004) 
which is based on s imilar principles as those applied by 
iBenson et al.l ||2002^. In their second mode l, they identi¬ 
fied sa tellites with the most massive subhalos. IZentner et alJ 
^200S^ found that in both cases, the satellite systems had a 
planar distribution similar to that in the Milky Way and 
argued that the degree of central concentration of the satel¬ 
lite systems plays an important role in this result. They 
also showed that the population of subhalos as a whole is 
anisotropic and preferentially aligned with the major axis of 
the triax ial halo. _ _ 

Like IZentner et al.l ^2^^!^^ , our study employs a semi- 
analytic model to follow the formation of the visible satel¬ 
lites. In this respect, bo t h the se studies are quite different 
from that of iKang et alJ (l200fil) who based their conclusions 
p urely on dark matter p articles. Our model differs from that 
of IZentner et al.l (l200!ih in several important respects. Our 
simulation codes and methods for identifying substructure 
are different. While they considered three halos specifically 
chosen to lie on a hlament, we used 6 simulations randomly 
chosen from a large cosmological volume. The biggest differ¬ 
ence, however, concerns the semi-analytic models used in the 
two studies. While both of them give a reasonable match to 
several observed properties of the Milky Way’s satellites, our 
semi-analytic model has been applied and teste d much more 
extensively than that of iKravtsov et al ] ll2004. The mode l 
we use is base d on the GALF D RM cod e of Gole et ^ (l200fi) 
as extended bv iBenson et al ] j2002bl) . This model has been 
shown to give an acceptable account of many properties of 
the galaxy population as a whole including the luminosity 
function in various passbands, from the UV to the far in¬ 
frared, and in various environments, distributions of colour, 
size and morphological type, etc. The model is also rela¬ 
tively successful at matching the properties of galaxies at 



Ntot 

(10®) 

Ahr 

(10®) 

-^vir 

kpc) 

A^vir 

(10®) 

ghl 

14.6 

12.9 

110 

1.07 

gh2 

18.1 

16.2 

131 

1.74 

gh3 

18.0 

16.2 

170 

3.73 

gh6 

25.5 

22.2 

169 

3.76 

gh7 

19.2 

17.3 

156 

2.99 

ghlO 

13.4 

12.1 

133 

1.86 


Table 1. Parameters for the six A-body halo simulations. For 
each halo (column 1), (2) shows the total number of particles in 
the simulation box, (3) the number of high resolution particles 
in the simulation, (4) the virial radius of the halo in h~^ kpc 
defined as the distance from the centre at which the mean interior 
density is 178pcrit i and (5) the number of particles inside the virial 
radius. In all cases, the simulation cube has co-moving length of 
35.325 h“'^Mpc, and the mass of the high resolution particles is 
2.64 X 

high-redshift, as discussed in iBa.iigfi et a,lJ (|2QQ^. Finally, 
the two studies use somewhat different methods to quantify 
the distribution of Milky Way satellites and to compare the 
results with the observations. On the whole, the conclusions 
of the two studies are consistent although there remain some 
differences as we discuss in Section 0| 

The remainder of this paper is organised as follows: in 
Section |21 we outline the methods used; in Section o we 
present our results which we interpret in Section r8.2l in Sec¬ 
tion 2] we discuss the implications of our findings. 

2 SIMULATIONS AND GALAXY 

FORMATION MODEL 

We have analyzed 6 high-resolution A-body simulations of 
galact ic-size dark matter ha los carried out with the GADGET 
code iSoringel et alJ^2^^1a^ . The halos, chosen to have a 
mass ^ IO^^Mq, were otherwise randomly selected from 
a large cosmological simulation of a cubic region of side 
35.325 h~^ in a flat ACDM universe (with flm = 0.3, h = 
0.7, (78 = 0.9). The simulation was executed a second time 
adding “high resolution” (i.e. small mass) particles, and ap¬ 
propriate small scale power in the initial conditions, to a re¬ 
gion surrounding the halo under consideration. These sim- 
ulations have been studied extensively in previous papers 

llPower et alJl2nn!i . llTa.vashi et al 11^0041 iNa.varro et al 11^00411 

and we refer the reader to those papers for specific details of 
how the simulations were carried out. Table 0 summarizes 
the important parameters of the simulations. 

We identified bound sub structures in the sim ulation 
using the algorithm SUBFINDjjSgnngeLe^L l2nnitJ) . First, 
“friends-of-friends” groups JPaU^etrdT^S^ are found by 
linking together particles whose separation is less than 
0.2 times the mean interparticle separation, correspond¬ 
ing roughly to particles within the virialized region of the 
halo. SUBFIND then identifies substructures within these ha¬ 
los based on an excursion set approach, using the spatial 
and velocity information for each particle in order to define 
self-bound objects. 

For each halo, we generate a complete merger history, 
identifying all progenitor and descendant halos, as described 
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in iHellv et, alJ (l2003f) . The semi-analytic galaxy formation 
model is calculated along each branch of the merger tree. 
This i s bas ed on the model descri bed in detail in ICole et alJ 
i200(tl and iBenson et al ] J2002ljl . The model includes the 
following physical processes: (i) the shock-heating and virial- 
ization of gas within the gravitational potential well of each 
halo; (ii) radiative cooling of gas onto a galactic disk; (iii) 
the formation of stars from the cooled gas; (iv) the effects of 
photoionization on the thermal state and cooling properties 
of the intergalactic medium; (v) reheating and expulsion of 
cooled gas through feedback processes associated with stella r 
winds and supernovae explosions ('see lRensnu et '^ l2003bll : 
(vi) the evolution of the stellar populations; (vii) the effects 
of dust absorption and radiation; (viii) the chemical evolu¬ 
tion of the stars and gas; (ix) galaxy mergers (which, de¬ 
pending on the violence of the merger, may be accompanied 

S irbursts and the formation of a bulge - see lBaugh et alJ 
; (x) the evolution of the size o f the disk and bul ge. 
Our m o del d iffers from that of ICole et alJ i200fll and 
iHellv et alJ (l20n.'t> in the way in which galaxy mergers are 
treated. In the current model, the positions of satellite galax¬ 
ies and the time when they merge is determined by using in¬ 
formation from SUBFIND. Central galaxies are placed on the 
most bound particle of the most massive subgroup in the 
halo. (SUBFIND identifies the background mass distribution 
of the halo as a separate subgroup, so this is generally a ro¬ 
bust way to define the centre of the halo.) Satellite galaxies 
are placed on the descendant subhalo of the progenitor halo 
in which they formed. If the subhalo ceases to be identified 
by SUBFIND at some later output time, we continue to trace 
its constituent particles and place the galaxy at the centre 
of mass of this group of particles. A galaxy is considered to 
have merged onto the central galaxy if its distance from the 
central galaxy is less than the spatial extent of the set of 
particles it is associated with. 

An overview of the results of our semianalytic model as 
regards the evolut ion of the gal a xy pop ulatio n as a whole 
may b e found in iBenson et al.l (l2003a^ and iBaugh et alJ 
i2004l while results r elevant to the satellit es of the Milky 
Way may be found in IBenson et alJ (|2002j. 


3 RESULTS 

We begin by quantifying the shapes of dark matter halos 
in the simulations and sub-sytems within them. We then 
interpret the results in terms of the formation histories of 
the halos and their subsystems. 


3.1 The morphology of halos and their subsystems 

The semianalytic model applied to the N-body simulations 
provides the position and internal properties of the central 
galaxy in each halo and its satellites. According to the semi¬ 
analytic model, three of the central galaxies are spirals and 
three e llipticals. For the pur pose of comparing with the anal¬ 
ysis of lKroupa et alJ (|200^, we select the 11 most massive 
satellites in each halo within a distance of 250 kpc from the 
central galaxy. We calculate the moment of inertia tensor of 
this satellite subsample, weighting each object equally, and 
obtain the principal axes of the distributions. 

Fig. m shows three orthogonal projections along the 
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Figure 1. Projections of the positions of the 11 most massive 
satellites within 250 kpc of the central galaxy along the principal 
axes of the inertia tensor in simulations gh6, gh2 and gh3. 


principal axes of the satellite systems in three of our six sim¬ 
ulations. The other three are very similar. The figure reveals, 
remarkably, that the loci of the 11 most massive satellites 
define a thin, disk-like structure around the central galaxy. 
As we show below, in most cases, the satellite structure is 
aligned with the major axis of its triaxial host dark matter 
halo. 

The distribution of the visible satellites differs signif¬ 
icantly from the distribution of the dark matter substruc¬ 
tures identified by SUBFIND. Fig.|^is analogous to Fig.^but 
the points plotted now correspond to the most massive 200 
substructures found within 250 kpc of the central galaxy. 
The projections are along the principal axes of the inertia 
tensor of the substructure systems. It is evident that the 
distribution of substructures is much less anisotropic than 
that of the satellites in Fig. 0 

The eigenvalues of the diagonalized inertia tensor are 
proportional to the rms deviation of the x, y and a co¬ 
ordinates relative to the principal axes. Denoting the ma¬ 
jor, intermediate and minor axes by a, b and c respectively 
{a > b > c), the flattening of the system may be quantified 
by the ratio s c/a a nd b/a. The early N-body simulations of 
iFrenk et alJ dlQSfjl showed that CDM halos are triaxial and 
recent work in dicates that c/a= 0.7 ± 0.17, and b/a> 0.7 

^Buiio^lio^. 

The axial ratios, found by diagonalizing the moment 
of inertia tensor, of the dark matter halos and various sub¬ 
systems of objects within them are plotted Fig. H Fig. EK 
shows that our simulated halos have axial ratios consistent 
with those found in previous simulations and tend to con¬ 
gregate near the top right of the panel corresponding to 
nearly spherical objects. With the exception of an outlier, 
gh6, which is significantly prolate, this is also the region 
populated by the massive subhalos. 

The axial ratios of the systems consisting of the 11 most 
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Figure 2. Projections of the positions of the 200 most massive 
dark matter substructures within 250 kpc of the central galaxy 
along the principal axes of the inertia tensor in simulations gh6, 
gh2, and gh3. 


massive visible satellites are plotted in Fig. int- The triangles 
correspond to our full semianalytic model and the squares 
to a variant in which the early reionization of the intergalac- 
tic medium is not included. The satellite systems in the two 
models have similar flattening because more then 80 % of 
the subhalos that host the brightest satellit es in the two 
cases a re the same. However, as discussed bv iBenson et alj 
Jiooii), neglecting the effects of reionization leads to an 
overprediction of the number of faint galaxies, including 
satellites in the Milky Way. Whether reionization is included 
or not, the satellites in our simulations cluster around the 
location of the Milky Way data marked by a cross in Fig. Et. 
This is the main result of our analysis: the flattening of the 
satellite system in our simulations is in excellent agreement 
with that of the Milky Way satellite system. 

It is clear from Fig. EK and Fig. Et that the most mas¬ 
sive visible satellites inhabit a biased subset of subhalos. To 
explore the origin of this bias, we select two subsets of sub¬ 
halos: the 11 most massive subhalos at 2 = 0 and the 11 
subhalos which had the most massive progenitors prior to 
being incorporated within the virial radius of the main halo. 
The flattening of these two systems is compared in Fig. Et- 
The figure shows that the crucial factor in establishing a 
highly flattened system is not the final mass of the subhalo 
but the mass of the largest progenitor. It is the latter that 
correlates well with the final stellar mass or luminosity of 
the visible satellite, as shown in Fig. BK Here we plot the 
stellar mass of each satellite galaxy against the mass of its 
largest progenitor. This strong correlation is a result if the 
GALFORM model readily making the most luminous galaxies 
in the most massive progenitro halos. In contrast, Fig. EJ) 
shows there is no correlation between the stellar mass of each 
satellite and the mass of its host substructure. This is due 



Figure 3. Minor to major (c/a) versus intermediate to major 
{h/a) axial ratios. Since a > b > c, the upper left triangular 
half of this plot cannot contain any points. Along the diagonal lie 
prolate objects and along the right vertical axis lie oblate objects. 
The numbers inside each symbol identify the simulated halo. The 
’X’ i ndicates the axial ra tios of the Milky Way’s satellite system 
from iKronna et, a] .11200,4 Only data out to a radius of 250 kpc is 
used in all cases. Panel (a) compares the axial ratios of the dark 
matter halos (triangles) with those of the system consisting of 
the dark matter substructures (squares). Panel (b) compares the 
axial ratios of the systems made up of the 11 most massive visible 
galaxies in models with and without early reionization (triangles 
and squares respectively). Panel (c) compares the axial ratios 
of the systems of 11 most massive substructures with those of 
the systems consisting of the 11 substructures that had the most 
massive progenitors. Panel (d) compares the axial ratios of the 
systems made up of the 11 most massive visible satellites in our 
full model (triangles) with those of the same systems but with the 
radial distances of each satellite normalized to a common value. 
Thi s same operation is peformed on the Milky Way data and so 
the lKroupa et ^l2Q0^ point moves slightly 


to the subhalos having been subjected to various amounts 
of tidal stripping. 

Comparison of Fig. Et and Fig. Et indicates that the 
flattening of the systems consisting of the 11 most massive 
visible satellites and the 11 subhalos that had the most mas¬ 
sive progenitors are very similar. This is an important result 
because it demonstrates th at our main c o nclus ion regard¬ 
ing the compatibility of the lKroupa et ahl ^200!^^ data with 
the CDM cosmology does not depend on the details of our 
semianalytic modelling of galaxy formation. So long as the 
brightest satellites form in those subhalos with the most 
massive progenitors, our conclusions stand. 

With only 11 satellites in our main samples, the possi¬ 
bility that the our estimate of the moments of inertia might 
be biased by the presence of outliers is a concern. We inves¬ 
tigate the sensitivity of our results to outliers by scaling all 
radial positions to a common value while keeping the angles 
of each radius vector fixed. The axial ratios of the rescaled 
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Least massive galaxies 



^^@lo(^prog halo/^ ^ ^o) subhalo/^ ^0) 


Figure 4. (a) Stellar mass in each satellite galaxy as a function 
of the mass that its largest progenitor halo had before becoming 
incorporated into the main halo, (b) Stellar mass in each satellite 
galaxy as a function of the mass of the substructure the galaxy 
currently resides in. The squares show results for the 11 most 
massive satellites in all 6 simulations, while the dots show results 
for the less massive satellites. The overlap in panel (a) is largely 
due to the smallest satellite mass varying between simulations 

data are compared to the axial ratios of the 11 most massive 
satellites in Fig. Hi. Rescaling the satellite radial distances 
scatters our estimates of axial ratios somewhat but does not, 
on average, lower the overall flattening of the systems. As 
shown in the figure, rescaling the Milky Way data in the 
same way also has a small effect on the axial ratios. 

Finally, we consider the connection between the highly 
anisotropic distribution of satellites and the orientation of 
their host dark matter halo. Consider the vector pointing 
along the major axis of the distribution (i.e. along Osat). 
Let 9 denote the angle between this vector and a vector 
pointing along the major axis of the halo, udm. For our 
six simulations, we find that cos(^) equals to 0.768, 0.979, 
0.702, 0.747, 0.387, 0.942 for galaxy halos ghl, gh2, gh3, gh6, 
gh7 and ghlO respectively. Thus, apart from gh7, there is a 
very strong alignment between the major axis of the disk¬ 
like satellite systems and the major axis of the parent dark 
matter halo. In the Milky Way, the major axis of the satel¬ 
lite disk-like structure is perpendicular to the galactic disk. 
Thus, if our galaxy resides in a dark matter halo similar to 
those that we have simulated, then the disk must be aligned 
such that its normal vector points in the direction of the halo 
major axis. Interestingly, this is exactly the alignment ob¬ 
served in recent gas dynamical simulations of the f ormat ion 
of spiral galaxies bv iNavarro. Abadi fc Steinmetd J2nn4h . 

3.2 Interpretation 

The highly anisotropic distribution of satellite galaxies in 
Milky Way type systems is a somewhat surprising outcome 
of galaxy formation in a CDM universe. This is particularly 
so in view of the fact that the population of subhalos as a 
whole is much less anisotropic and has axial ratios similar 
to those of the halo dark matter. The key to understand¬ 
ing the origin of the anisotropic satellite distribution lies in 
the connection between halos and the cosmic web and, in 
particular, in the way in which satellites are accreted onto 



Figure 5. The formation of a galactic halo and its satellites. The 
points show a random 1% of the dark matter particles that end up 
in the main halo and the red circles the positions of the 11 most 
massive satellites that end within 250 kpc of the main galaxy by 
the present day. The scale of each plot is indicated by the red 
line which has a co-moving length of 400 kpc. The initial collapse 
produces a 2D structure — a large pancake of dark matter. 





















6 N. /. Liheskind et al. 


the main halo. Fig. illustrates the anisotropic nature of 
satellite accretion. The dots show a random 1% of the dark 
matter particles that end up in the main halo at the final 
time. The red circles mark the locations of the most massive 
progenitors of the 11 satellites that have the largest stellar 
mass at the final time. Rather than originating isotropically, 
those halos destined to become bright visible satellites are 
accreted primarily along one or two of the cosmic web fila¬ 
ments. 

The first panel in the figure illustrates the highly 
anisotropic collapse typical of CDM structures on galac¬ 
tic scales. Careful inspection of the time evolution of the 
system shows that the collapse occurs first along 2D sheet¬ 
like structures which subsequently wrap up into filamentary 
streams of dark matter. By 2 ~ 4.2, these filamentary “high¬ 
ways” along which proto-satellite galaxies form are well es¬ 
tablished. The filaments are generally thicker than the locus 
of the largest proto-galactic halos which tend to concentrate 
towards the central, densest parts of the filament in a near 
1-dimensional configuration. As the most massive halo pro¬ 
genitors collapse to form the main galaxy, this alignment is 
largely preserved. Smaller halos are more widely scattered 
across th 
strength 

addition, they are often accreted over a longer period and 
from a larger range of directions. Their distribution, now 
lacking a preferred orientation, ends up being much less 
anisotropic than that of the most massive halos. Whether 
reionization is included or not, satellite galaxies in the semi- 
analytic model form preferentially in the subhalos with the 
most massive progenitors and thus inherit their highly flat¬ 
tened configuration. 


e thick filaments, reflecting their weaker cluste ring 

). In 


4 DISCUSSION AND CONCLUSIONS 

We have shown that the, at first sight surprising, flattened 
distribution of satellites in the Milky Way is the natural out¬ 
come of the anisotropic accretion of matter along a small 
number of filame nts, characteri s tic of halo formation in the 
CDM cosmology. iKrnnna et ail (l20n!jl reached the opposite 
conclusion, that the observed satellite distribution is incom¬ 
patible with the CDM model, because they neglected the 
fact that the satellites do not trace the distribution of halo 
dark matter but form instead in the most massive halos 
whose spatial distribution is biased. 

Our re s ults are not directly comparable to those of 
iKang et al.l (1200 pfl who also attempted to interpret the flat¬ 
tened distribution of Milky Way sa tellites with the ai d of 
high-resolution N-body simulations. iKang et alJ (l200ljl as¬ 
sumed that the satellites follow the dark matter distribution 
in the halo an d did not consider th e formation sites of satel¬ 
lites in detail. IZentner et alJ ||200^, on the other hand, im¬ 
plemented a semi-analytic model similar to our s in N-body 
simulations also similar to ours and those of iKang et alJ 

m . 

Our results ar e broadly consistent with those of 
IZentner et al.l i2nn,^ . Unlike them, we did not choose halos 
specifically lying along filaments but selected them at ran¬ 
dom from a large cosmological simulation. In the event, three 
of our halos would, according to our semianalytic model, 
host spiral galaxies and the other three elliptical galaxies. 


One d ifference between the two studies is that IZentner et ^ 
^200S^ found an acceptable match to the Milky Way satellite 
distribution both in their semianalytic model and in a model 
in which the satellites are identified with the most massive 
subhalos at the final time. We have shown that the distri¬ 
bution of the latter is not as flattened as the distribution 
of Milky Way satellites. The crucial factor is not the final 
mass of the halo which is affected by tidal stripping, but the 
mass of the largest progenitor before it is accreted into the 
main halo. Indeed, if the satellites are identified with the 
halos that had the largest progenitors, then their flattened 
distribution is very similar to that of the satellites identified 
by our semianalytic model. Thus, our conclusions are inde¬ 
pendent of the details of our galaxy formation modelling. 

As was also found by IZentner et al.l (l200?ji . the ma¬ 
jor axis of the flattened satellite distribution in our sim¬ 
ulations points close to the direction of the major axis of 
the parent halo. This alignment reflects the preferential 
accretion of mass onto the halo along the dominant fila¬ 
ment. An important consequence of this result is that if the 
Milky Way resembles the systems we have simulated, then 
the Galactic disk should lie in the plane perpendicular to 
the major axis of the halo because the observed satellite 
system itself is perpendicular to the Galactic disk. Inter¬ 
estingly, this is exactly the configuration found in the re- 
cent gasdynamical simulations of sp iral galaxy formation by 
iNavarro. Abadi fc Steinmet d ll2004ll . 

The satellite alignment that we have found in our simu- 
lations is almost certainly related to the “Holmberg effect,” 
iHolmberdllQbflll the observation that the satellites of exter¬ 
nal galaxies within a projected radius of Vp ~ 50 kpc tend 
to lie preferentially in a cone along the galaxies’ minor axis, 
avoiding the equatorial regions. To test this observation re¬ 
quires a larger number of simulations than those we have 
performed. Similarly, our current simulations are inadequate 
to test the extension of the Holmberg effect uncovered by 
IZaritskv et alJ (|^93) from a study of isolated spirals which 
also revealed an excess of satellites along the minor axis of 
the galaxy, now out to project ed distances of 50 0 kpc. 
A similar result was found by ISales fc Lambaa (l2004^ from 
a much larger sample of galaxies drawn from the 2 degree 
field galaxy redshift survey. They too found an anisotropic 
distribution for Vp < 500 kpc, but only for satellites moving 
with a velo city relative to th eir host of An < 160 km s 
In contrast, iBrainer^ ll2004^ found the opposite effect in a 
sample of satellites from the Sloan Digital Sky Survey, an 
alignment along the major axis at small radii (rp < 100 kpc) 
and an isotropic distribution beyond. 

Although our simulations are not large enough to study 
the distribution of satellites beyond the inner 250 kpc of the 
galactic centre, it seems likely that the anisotropic distribu¬ 
tion of satellites will continue out to larger separations. We 
intend to study this problem in a larger set of simulations. 

In summary, we have found that flatte ned distribu¬ 
tion o f the Milky Way satellites, first noted byjLgidei^Bel 
iW^), and most recently highlighted by iKrouna et al 
I 200^, turns out to have a simple explanation in the con¬ 
text of structure formation in the GDM model. It is merely 
a reflection of the intimate connection between galactic dark 
matter halos and the cosmic web. 
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